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Abstract. Mode entanglement exists naturally between regions of space 
in ultra-cold atomic gases. It has, however, been debated whether this 
type of entanglement is useful for quantum protocols. This is due to a 
particle number superselection rule that restricts the operations that can 
be performed on the modes. In this paper, we show how to exploit the 
^ mode entanglement of just a single particle for the teleportation of an 

Q unknown quantum state of a spatial mode. We detail how to overcome 

the superselection rule to create any initial quantum state and how to 
perform Bell state analysis on two of the modes. We show that two of 
the four Bell states can always be reliably distinguished, while the other 
two have to be grouped together due to an unsatisfied phase matching 
__ condition. The teleportation of an unknown state of a quantum mode 

thus only succeeds half of the time. 

! 1 1 Key words: Mode entanglement, quantum teleportation, superselec- 

C2 tion rule 

, ^ 1 Introduction 

t-H Entanglement is a key resource in many practical applications using quantum 

mechanics pp. Usually entanglement is thought to exist between the degrees 
of freedom of two, or more, well localised quanta, such as photons or massive 

j> particles. When particles are well separated, they are distinguishable from one 

£f) another and can thus be assigned labels. Entanglement between the particles 

is then well defined since their Hilbert space has a tensor product structure. 
However, if the de Broglie wavelengths of identical particles begin to overlap, 
the particles become indistinguishable and one can no longer assign a label to 
each particle. The concept of standard particle entanglement breaks down as 
the Hilbert space no longer has the required tensor product structure; it is a 
projection onto the symmetric or antisymmetric subspaces (depending on the 
particle statistics) [2]. 

While entanglement between indistinguishable particles can still be correctly 
defined by taking a set of detectors into account [3], another method to recover 
a tensor product Hilbert space is to move into the formalism of second quantisa- 
tion. Here, in the so called Fock basis, one defines a complete set of single particle 
states and counts the number of excitations in each. The corresponding mode 
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structure has a tensor product Hilbert space and hence entanglement of modes is 
a meaningful concept. Entanglement can therefore exist between modes occupied 
by particles [214] . The modes can be energy eigenmodes or perhaps more rele- 
vant to quantum communication or information processing protocols are spatial 
modes. In the following section, we will discuss in more detail the mode structure 
for a simple many-body system. The particles may be massless, i.e. photons, or 
massive, such as those found in ultra-cold gases. Mode entanglement of photons 
has been considered in a number of works |5l(i| and is usually limited to the 
single photon regime [7 8 9 10J. The experimental confirmation of the mode en- 
tanglement of a single photon via a Bell like test was obtained in 2004 ,TT] and 
the multipartite entanglement of one photon distributed between four optical 
modes was characterized using uncertainty relations in an interferometric setup 
in 2009 [12]. 

For massive particles, the existence of entanglement between spatial modes 
becomes less clear [13] . This is due to a particle number superselection rule 
[14115] , that forbids isolated systems from existing in a superposition of eigen- 
states of different mass. Hence, for modes occupied by massive particles the 
system density operator, p, cannot contain any off-diagonal terms that connect 
states of different particle number. Any measurements made on the spatial modes 
are also restricted to the subspace of fixed particle number. Since the correla- 
tions of entanglement are locally basis independent, they can be confirmed via 
for instance a Bell inequality, where one should measure each subsystem locally 
in at least two basis. For spatial modes, the particle number is one such basis, 
however to measure in a second, rotated basis, measurements of superpositions 
of different numbers of particles are required. For an isolated system, such su- 
perposition measurements are forbidden due to the number of particles in the 
system being fixed. Until recently it was unresolved how to measure modes of a 
massive bosonic field in any way other than in the particle number basis. And 
this is why it has been debated [16| whether mode entanglement is as 'genuine' 
as particle entanglement or whether it is just a mathematical feature of the 
quantum state. 

Recent research has, however, shown that superselection rules [17 118119 ] can 
be overcome locally by using a suitable reference frame |17l20j . The reference 
frame required to rotate spatial modes away from the particle number basis is 
a coherent reservoir of particles such as a Bose-Einstein condensate (BEC). By 
using such a reservoir, one can create, at least in principle, superpositions of dif- 
ferent numbers of particles [21120] . Thus, it is predicted that mode entanglement 
of a single massive particle is, indeed, as genuine as particle entanglement |22| . 
Specific schemes for Bell inequality tests of mode entanglement of one [23 25124] . 
and more [26] . massive particles have been given. Another question is whether 
mode entanglement can be used for quantum communication and quantum in- 
formation processing. This is an interesting point, as mode entanglement occurs 
naturally in coherent ultra-cold bosonic gases [27 28 29 30 31 . If one could har- 
ness this entanglement for practical applications, it could by-pass the need to 
create complicated entangled states manually. 
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A first step to understanding the usefulness of the mode entanglement of 
just a single massive particle came in a recent paper |32j . which gave a scheme 
for implementing the quantum dense coding protocol. Dense coding |33j allows 
the transmission of two classical bits of information via one qubit; in order to 
achieve the full quantum channel capacity a maximally entangled Bell state is 
initially required and full Bell state analysis is needed. It was shown [35] that 
the linear photonic dense coding channel capacity could be achieved without a 
BEC reservoir and that with a BEC reservoir the full quantum dense coding 
protocol could be implemented. 

While it was briefly mentioned in [35] that the teleportation of a quantum 
state of a spatial mode should also be possible, no detailed scheme was given. 
It is the aim of this paper, to provide such a scheme. Quantum teleportation 
[34 is a key protocol in quantum information science and has been fruitfully 
demonstrated in experiments with a number of different physical implementa- 
tions 35 36 37]. In [35J, a qubit (or optical mode) consisting of the vacuum and 
one photon states was teleported using the entanglement of a single photon in a 
superposition of two spatial modes. Even though the creation of a superposition 
of different numbers of photons is not strictly forbidden by a superselection rule 
as it is with massive particles, there is still the problem of how to keep track 
of the phase between the two different photon number states. Lombardi et al 
[35] solved this by teleporting a mode that was entangled to another one. In 
other words, the single photon was coherently distributed across the mode to be 
teleported and an ancilla mode. This ancilla optical mode actually played the 
role of a phase reference in a similar way to how the BEC reservoir will play the 
role of a phase reference in our scheme here. 

One motivation to study a teleportation scheme using the mode entangle- 
ment of a single massive particle is to allow for direct comparisons of quantum 
information processing with different physical systems. Moreover, from a funda- 
mental viewpoint, it is interesting to clearly demonstrate that mode entangle- 
ment of massive particles is, in principle, useful entanglement - and can actually 
be used in a very similar way to particle entanglement - for quantum information 
processing despite the particle number superselection rule. In particular, in this 
paper we will see that in contrast to the dense coding scheme of [35] , where the 
full channel capacity was attainable, here an arbitrary state of a quantum mode 
is only reliably teleported half of the time. This is due to an additional phase 
locking criterion that arises from having a total of three modes in the telepor- 
tation scheme as opposed to the two modes that are required for dense coding. 
This illustrates subtle intricacies of teleportation with mode entanglement of a 
massive particle that are not present in the particle entanglement case. 

We begin in the subsequent section, by reviewing the concept of mode entan- 
glement. In section (13]), we detail how to teleport the unknown state of a spatial 
mode using a single massive particle distributed coherently across two spatial 
modes. We begin next by introducing the concept of spatial mode entanglement 
in more detail. 
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2 Entanglement of spatial modes 

We will start by detailing the mode structure of a bosonic system and then use 
a simple example of mode entanglement to illustrate the differences between it 
and particle entanglement. Finally, we will explain some results concerning mode 
entanglement of Bose gases at zero and finite temperatures. Note that we will 
consider in this paper only mode entanglement of bosonic fields; for a discussion 
of mode entanglement in fermionic systems see, for example, this paper |39) by 
Aharonov and Vaidman. 

Consider a confining volume, V, whose energy eigenmodes, labelled by fc, can 
be 'excited' by applying creation operators, a k , on the vacuum state, a\\vac) = 
|0i02...1fe0fc + i...), where [a,t,a|] = 5 k i- An excitation is 'a particle' with corre- 
sponding energy E/~ = fouk, where cok is the frequency of the fc-th energy eigen- 
mode. As mentioned already, instead of describing the system using its energy 
modes, it is often desirable to use a different set of modes, such as the spatial 
modes. The description of the system in space is obtained by a transformation 
via the energy eigenf unctions, a\\vac) = J dx <fik (x) (x) \vac), where <fik(x) is 
the fc-th energy eigenfunction and ip'(x) creates a particle at point, x, in space. 
It follows that fy(x), i)(x')] = S(x — x'). Populating an energy mode with a par- 
ticle is thus equivalent to populating all spatial modes, i.e. all points in space, 
in a superposed manner. 

To illustrate the differences between particle and mode entanglement, we 
take two non-interacting bosons trapped in a confining volume at zero temper- 
ature. In first quantisation, i.e. in the language that one would use to describe 
particle entanglement, the wavefunction is the symmetrized product, ^12(2:, y) = 
■^={^\{x)(j)2{y) + 4>\{y)(l)2{x)), where (f>(x) is the ground state of the confining 
potential. No entanglement exists between the particles, since indistinguishabil- 
ity forbids us from assigning to any particle a specific set of degrees of freedom. 
In other words, the state space of the two particles is a projection onto the 
symmetric subspace, whereas a tensor product Hilbert space, H = Hi <§5 7^2, is 
required to define entanglement between the subsystems. 

Conversely, in second quantisation one can define a pair of spatial modes, A 
and B, where each mode occupies half the confining geometry. Since both the 
particles are coherently distributed over these modes, the system is described by 
the entangled state, 




\vac) 



2 




i(|20) + V2|ll> + |02» 



(1) 



where 'ip x = f x dxg(x)fo(x) creates a particle in mode X = A, B (g(x) is a 
so called detector profile that gives weighting to the points in space), |mn) = 
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\rn)A ® I"-) b span the state space H = Ha ®Hb and m denotes the number of 
particles in mode A and n the number particles in mode B (with m + n = 2). 
From this example it is clear that entanglement is contingent on the choice of 
modes [JT], but provided that like here a suitable choice is made, investigat- 
ing entanglement between distinguishable modes circumvents the difficulties of 
defining entanglement between indistinguishable particles |40j . 




Spatial mode entanglement of an non-interacting BEC at zero temperature 
was first considered by Simon [27] . He and others [H] found that mode entangle- 
ment existed between regions of space if the gas had an uncertainty in particle 
number below a given level. That is, if the gas is best described by a coherent 
state, | a) (or mixtures there of), there is no entanglement between spatial modes. 
On the other hand, if a uniform gas is of a fixed particle number, N, there is 
\ log 2 N amount of entanglement (as measured by the von Neumann entropy), 
between two equal sized modes. Anders et al [2 8) used a thermodynamical en- 
tanglement witness to show that spatial mode entanglement only exists across 
an entire Bose gas below the critical temperature for BEC. Spatial entanglement 
between two and more modes of an interacting Bose gas at finite-temperature 
with a fixed (but possibly unknown) number of particles was considered by Goold 
et al. |31j . who demonstrated that mode entanglement is present in a gas when 
the coherence length of the particles extends over the modes. More specifically, 
there is a direct link between single-particle reduced density matrix [43] (i.e. the 
visibility of interference fringes) between different regions of the gas and spatial 
mode entanglement. This means that the natural mode entanglement of a BEC 
has already been detected in experiments such as [U], albeit indirectly. More 
recently, mode entanglement generated by a single exciton was predicted to exist 
between different sites in the photosynthetic FMO complex [45], (as a result of 
the experimentally verified quantum coherence in the molecule), demonstrating 
that this type of entanglement is relevant even in biological systems. 

3 Teleportation of a quantum state of a spatial mode 
using a single massive particle 

Since mode entanglement exists naturally within many systems, it is important 
to ask whether, at least in principle, this type of entanglement is useful entangle- 
ment. We will address the question by providing a scheme to show that a single 
massive particle that is coherently distributed over two spatial modes can, at 
least theoretically, be used to teleport the unknown quantum state of a spatial 
mode perfectly half of the time in spite of the superselection rule. 

In the following, we will first introduce a Hamiltonian whose parameters can 
be switched on and off to perform the gates. We will show that in order to 
overcome the superselection rule and to rotate the modes away from the particle 
number basis, a Bose-Einstein condensate should be used as a particle reservoir 
and also as a phase reference throughout the teleportation protocol. We will 
end by introducing the quantum circuit that allows to perform the teleportation 
protocol. 
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The standard qubit teleportation protocol [34] between two parties, A for 
Alice and B for Bob, can be split up into four parts: 

(i) preparation stage: the preparation of an unknown quantum state by a 
third party, Charlie, and also the distribution of an entangled Bell state between 
Alice and Bob, 

(ii) Bell state measurement by Alice: the prepared qubit is passed to Alice 
who then makes a Bell state measurement on this and her portion of the entan- 
gled state and records the measurement outcome, 

(hi) transmission of classical information: Alice sends Bob two bits of clas- 
sical information that indicate which of the four Bell states was measured, and 
(iv) single mode rotation by Bob: Bob uses the classical information to select 
which operation to perform on his qubit leaving him with the original unknown 
quantum state. 

We will now detail how to perform the four steps using mode entanglement of a 
single massive particle. 

The system consists of three spatial modes, a, A and B. Mode a will be placed 
in an unknown state, which will be teleported to mode B. Modes A and B will 
be in the maximally entangled state formed from a single particle. The system 
is described by the Bose-Hubbard model with additional coupling of each mode 
to the BEC reservoir. The Hamiltonian is written as follows 

Hbh = -^-(4>a4>b + V'bV'a) - -y^I^a + $a4>o) 

+ J2 Uini(ni-1)+ E **i- y($^re S +$ e >iX2) 

i=a.A,B i=a,A,B i=a,A,B 

The first two terms represent the coupling between the modes, A and B, and 
also modes, a and A. Modes a and B do not interact throughout the entire 
protocol. The coupling between the modes can be turned on and off by varying 
the tunneling matrix element J a A {Jab) by increasing/decreasing the height of 
the potential barriers between the wells by altering the intensity of the trapping 
laser. 

The parameter Ui is the onsite interaction term, which is continuously set to 
be much larger than the other energy scale in the Hamiltonian. This is achieved 
by ensuring that each potential is tightly confined in all three directions, as is 
the case for atomic quantum dots [21] . The resulting large nonlinear repulsive 
interaction between the particles means that only zero or one particles can exist 
in the mode at any instance. In other words, by maintaining a high repulsivity 
between the particles the modes are forced to behave like qubits. 

The next three terms in the Hamiltonian are the free energies Ei of the 
individual modes. The standard setting throughout the protocol will be E a = 
Ea — Eb, but to change the phases of the modes relative to one another, a 
potential bias can be applied to a mode using, for instance, a dispersive laser 
pulse. 
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The final set of terms in the Hamiltonian correspond to the coupling of each 
individual mode to the BEC reservoir. We consider a Raman laser set up, which 
couples the two different trapping states of the system and reservoir atoms [?Bj . 
The parameter, Q — J G?x</>o(x)<J r o(2 ; )^j is the effective Rabi frequency, where 
&q(x) is the wavefunction of the BEC, (f>o(x) is the wavefunction of the atom in 
mode a and J? is the usual Rabi frequency. 

3.1 Single-mode and two-mode gates 

We will now discuss the implementation of the single-mode and two-mode gates 
that will be used in the teleportation protocol. 

Single-mode phase gate: A phase gate, for instance a Pauli Z operation, can 
be applied to mode j by altering the energy, Ej, relative to the free energies 
for the other modes for a given time, t. The corresponding unitary operation on 
the mode is U = e l£j ™ j t , which if applied for t = ir/Ei performs the Z gate, 
|0) — > |0), |1) — > — During this operation, the couplings between the modes 
themselves and between the modes and the BEC are switched off. 



Single-mode number rotation gate: We would also like to rotate the state 
of an individual mode between the |0) and |1) eigenstates. To do this we need to 
couple to the BEC reservoir by applying the Raman pulse for a time, t. We take 
here a uniform non-interacting BEC described by a mixture of coherent states, 
Pbec — J Q 27r |f ||a|e l8 )(|a|e ie |, wnere M 2 = " is the average number of particles 
in the condensate and 9 is the condensate phase. For a nummerical analysis of 
a single mode, i.e. an atomic quantum dot, coupled to a superfluid reservoir of 
interacting particles see [37] , 

Note that when calculating it suffices to use the coherent state, ||a|e ), 
with one realisation of the phase 9 and to apply the twirling operator, T[pl] = 
/„" f^[/3i], to obtain the quantum state of laboratory as seen by Bob (here pl 
is the density operator for the total system and laboratory |20|17j ). Note that, 
unlike the dense coding scheme [32], Alice and Bob do not necessarily need to 
use the same BEC for their operations in this teleportation scheme. The unitary 
evolution governed by Hint = ^(^JVVes + ^J es V'i): transforms the occupation 
number of the modes as follows 

|0>- cob(^)IO) -*sin(^K«|l) 

|l>^c OS (^)|l>-isin(^)e--|0>. (3) 

Here we have traced out the BEC, which remains separable to the state of the 
modes as we assume it to have a high mean number of particles, n >> 1. By 
coupling the modes to the BEC for various times, t, single mode number rotation 
gates can be performed. 
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Two-mode c-phase type gate: We also need a two mode entangling gate 
to allow for Bell state analysis. A c-phase type gate between two of the modes, 
j =/= k G a, A, B, is performed as follows. Since we consider a large onsite inter- 
action strength throughout this protocol, the bosons behave like spin-polarised 
fermions. One can see this by associating the two particle number states, |0) and 
|1), of each mode with the up/down spin half degree of freedom. The bosonic 
Hamiltonian, H = —^(ipjipk + 4>k t l , j)' becomes equivalent to the quantum XX 
spin model, which in turn can be mapped via the Jordan- Wigner transforma- 
tion to H = —J(djCk + c\.Cj), where c x and cx are Fermionic creation and 
annihilation operators for mode, X — j, k that anti-commute, {cx, cy}+ = 5xy 
[49] . Once the system is in this regime, the barrier between the modes, j and 
k, is lowered for time, t — 7r/(2J), so that the particles exchange position. This 
ensures that the |11) term picks up a minus sign |11) — > — 111), due to the anti- 
commutation relations [50]. The other three states transform as, 1 00) —> |00) 
1 10) |01) and |01) -> 1 10). 



3.2 Circuit diagram for teleportation of a quantum state of a mode 

We now present the circuit diagram for the teleportation protocol, see Fig. ([IJ, 
referring to the four steps outlined at the start of this section. The three modes 
and the BEC are represented in the figure by the four horizontal lines in the 
diagram and the gates are represented as boxes. 

The preparation stage: a third party, Charlie, prepares a spatial mode (a) 



Preparation 



Bell state measurement 



B |0> 
A|l>- 

« |0>- 



U(,/2) 



BEC 



amy 



U(„/2) 



U*00 



U(,/2) 



IV2) 



Fig. 1. The circuit for teleportation of a quantum state of a spatial mode. The 
three modes and the BEC are represented in the figure by the four horizontal 
lines in the diagram and the gates are represented as boxes. The first portion of 
the circuit is the preparation of the entanglement via a single particle between 
Alice and Bob and the creation of an unknown state of the spatial mode a. The 
Bell state measurement on modes a and A by Alice is outlined in the centre 
portion of the circuit. The last part of the circuit details the final rotation on 
mode B by Bob. Further details will be discussed in the text. 



in an unknown quantum state by performing a single-mode number rotation 
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(green in the diagram [section (3.1 )]) and a single-mode phase gate (red in the 
diagram [section ( |3.1| ]) for the chosen amounts of time. Initially, mode a, is 
in the vacuum state \a) — |0) a , so that the two gates rotate it to the state 
| a) = a|0)+/3e J ( e+ ^|l), where a = cos (00, /3 = -£sin(0O and & = V^/Et/2 and 
4> = Eit. The potential barrier between modes a and A remains high through- 
out this stage. An entangled Bell state shared between Alice and Bob is also 
required for the protocol. To obtain such a state a single massive particle needs 
to be coherently distributed over the two modes - if the single neutral atom 
initially starts in mode A, lowering the barrier to mode B for time t = n/(2J) 
creates the state \^P + )ab = -^.(]10) + |01) ). 

The Bell state measurement: Spatial mode a is now passed to Alice, who 
performs Bell state analysis on modes, a and A. Bob is in possession of mode 
B. A measurement in the Bell state basis proceeds as follows: 

(1) Modes a and A should be within close proximity of one another and the 
potential barrier between them is high so that tunneling is fully suppressed. Mode 
A is rotated away from the particle number basis by coupling to a BEC reservoir 
via the Raman laser set-up for time t — n / (2fi^/n) as denoted by the yellow 



gate in the circuit diagram [see section (3.1 )]. Note that the BEC here does not 



necessarily need to be the same as the one used by Charlie to create the initial 
state of mode a. This results in the transformations, \Q)a ^(1^) — «e tff 1 1) ) a 
and \1)a — > ^(l 1 ) — is \0))a- For instance, the \ip+) a A Bell state transforms 

(without normalisation) to \01) aA + \10) aA -> |0) o (|l) - ie' i6 \0)) A + |l) o (|0) - 
ie \1))a- The other Bell states transform in a similar way. 



(2) Then, a c-phase type gate [see section (3.1 1] is applied by lowering the 
potential barrier between modes a and A for time t = n/(2J) so that the bosons 
swap positions. This is illustrated by the pale blue gate in the circuit diagram. By 
continuing the example in the previous step, one can see that this transformation 
works in the following way: |0)„(|1) - ie- w \0)) A + |1),(|0) -ie ie \l)) A -> \01) aA + 
\X0) aA - ie- w \00) aA + ^e ^e |ll) Q A = (|0) + ie w \l)) a (\l) - ie- ie \0)) A . 

(3) In order to rotate modes, a and A, to the particle number basis for the 
read-out, each should be coupled to the BEC reservoir for time, t — n/(2n^/K), 
as in step (2). For example, the state, (|0) + ie te \l)) a (\l} — ie~™\0))A, after c- 
phase type gate above becomes |00) a ^. Hence a measurement of zero particles 
in both modes means that the quantum state, \a) = a\0) + (3e l ( e+ ^\l) has been 
teleported to Bob. A measurement of |01) a A corresponds to the state, Z\a) = 
a|0) — f3e^ e+ ^\l), being teleported to Bob's mode. Note that the teleported 
states both also have phases which are correlated to the BEC, which means if 
they are to be used for other purposes then the same BEC would be needed to 
ensure phase matched conditions. 

However, if there is one particle in mode a after this final rotation, the tele- 
portation has not succeeded. This is due to the fact that after steps (l)-(3), the 
states, |0 ± ) a A, are rotated to |1) ((1 ± e l2e )|0) + (1 T e i29 )\l)) A /2. Since one 
does not know the phase, 9 of the condensate, one has to average over all re- 
alisations of it, i.e. apply the twirling operator to determine the quantum state 
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of the laboratory, which accounts for our ignorance. This leaves mode A in the 
maximally mixed state so that it is impossible to say which of the two operators, 
X and ZX, should be applied to recover the original state, \a), even in principle. 

Transmission of classical information and final single-mode rotation: 

Alice thus sends one of three messages to Bob according to the outcome of her 
measurement. Bob obtains the unknown quantum state half of the time upon 
the following transformations: 

\^ + )aA -> |00) aA -> tg, \i>-) aA -> |01) aA -> Z B 

\(f> ± )aA — > |l)(l|a <& Ia teleportation failed (4) 

where the first column corresponds to the Bell state Alice measured, the second 
column to the particle number measurements in mode a and A (and also the 
classical bits that she sends) and the third column to Bob's operation on mode 
B. The Z operation is applied to mode B by changing the phase of the mode [as 



in section (3.1 )] for time t = n/ Eb where Eb is the energy bias of mode B. This 
concludes the teleportation protocol for a unknown quantum state of a spatial 
mode. 



4 Discussion and conclusion 

In this paper, we have discussed the teleportation of an unknown state of a 
spatial mode using a mode entangled state formed from a single particle. We 
modeled the modes as qubits by considering tightly confined potentials, with 
zero or one particles representing the qubit degrees of freedom. We have given 
an explicit scheme for creating an arbitrary state of a spatial mode despite the 
particle number superselection rule, by coupling to a BEC reservoir. We have 
shown that one can reliably distinguish two of the four Bell states, with the other 
two grouped together. This is due to the fact that we do not know the phase 
of the BEC reservoir, which is still classically correlated to the state after the 
Bell state analysis. However, the teleportation of a general unknown (two-level) 
state of a spatial mode is still achieved half of the time, which would never be 
possible if we could not locally bypass the superselection rule in the first place. 

This contrasts with a dense coding scheme using the mode entanglement of 
a single particle (see [32]). In this scheme, one can always discriminate all four 
Bell states. The difference between the two schemes arises from the difference in 
symmetries of the encoding (or preparation) and decoding (or Bell state analysis) 
processes in the two protocols. In the dense coding protocol, there are only two 
modes and Alice can initially couple her mode to the BEC to rotate to one of 
the two states, l^) = ^(|00) ± e l20 |ll)). In doing so, the phase of these states 
becomes classically correlated to the phase of the BEC, 9. Alice then passes her 
mode to Bob, who performs Bell state analysis as was outlined above to find 
which Bell state Alice sent. Because the phase of the BEC is already present 
in j^) it cancels with the phase that is picked up again from the BEC during 
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the Bell state measurement. Thus all four Bell states are reliably distinguishable 
from one another. 

On the other hand, in the teleportation scheme, there are necessarily three 
modes. Initially the mode that will be subsequently teleported to Bob is rotated 
to an arbitrary state and becomes classically correlated to the phase of the BEC - 
this phase will stay present in the state of this mode as it is teleported to Bob and 
does not feature in the Bell state analysis. Bell state analysis is then performed 
on this and one mode of the entangled state. However, the phases that are picked 
up from the BEC in the Bell state analysis do not have any phases to cancel 
with, unlike in the case of the dense coding scheme. So while we have access to 
all of the gates that are required to perform the teleportation protocol, because 
we cannot know the phase of the BEC, even in principle, we do not teleport the 
spatial mode all of the time. This paper has illustrated how mode entanglement 
of massive particles does not behave exactly like particle entanglement despite 
the fact that one can locally overcome the particle number superselection rule. 
Future work should focus on extending this scheme to consider the teleportation 
of a general d > 2 level state of a spatial mode using the natural entanglement 
in BECs. 
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